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1. Introduction 

In the path integral formulation of fermionic field theories, one is forced to confront 
the determinant of the Dirac operator, while for spin-0 bosons one encounters the 
determinant of the gauged Laplacian. In the absence of gauge fields, and if the space- 
time geometry is not an issue, this determinant is an irrelevant constant. However, 
there is still a good deal of work to be done toward elucidating the dependence of 
these determinants on background gauge fields and gravitational fields. These deter- 
minants are related to the one-loop effective action via S'cff oc lndetZ?, and, in the 



'Present address [2009]: AIP, Washington DC; e-mail: sblau@aip.org. 
^Present address [2009]: School of Mathematics, Statistics, and Operations Research, 
Victoria University of Wellington, New Zealand; e-mail: matt.visser@msor.vuw.ac.nz; 
url: http: //homepages. msor rvuw.ac.nz/~visser/| 

t Present address [2009]: Theoretisch-Physikalisches- Institut, Friedrich-Schiller-Universitat Jena, 
e-mail: 07743 Jena, Germany, e-mail: wipf@tpi.uni-jena.de; url: http://www2.uni-jena.de/~p5anwi/ 



1 



June 16, 2009 6:34 WSPC/Guidelines-IJMPA analytic-4-arXiv 



2 S. K. Blau, M. Visser & A. Wipf 

guise of the one-loop effective action, have been the subject of considerable efforts 
dating back at least to the seminal works of Euler and Heisenberg, ^ Weisskopf, ™ 
and Schwinger. ^ 

Only in two dimensions is the situation reasonably well understood. The par- 
ticularly simple geometry of compact two-dimensional Riemann surfaces makes it 
possible to express the gauge-field-free Dirac determi nant in terms of Riemann 
theta functions and generalised Dedekind eta functions. I 4 | 5 | 6 | 7 | Q n ^ ne (-^her hand, 
in the plane the determinant for localized gauge fields is given by Schwinger's ^ 
famous result S e s oc / A^A^. Combining the effects of non-trivial topologies and 
non-zero field strength has been discussed in reference^. 

In higher dimensions, only particular background fields can be handled. For 
example, the effect of co nfor mal metric deformations on the effective action is dis- 
cussed in references and Schwinger El has considered the effective action 
for constant electromagnetic field strength and for a plane wave of electromag- 
netic radiation, both in flat four-dimensional Minkows ki space. Four-dimensional 
instanton solutions have also been considered. |12 | 13 | 14 [ 

In this paper we shall be interested in obtaining explicit analytic results. Ac- 
cordingly, we are forced to restrict our attention to particular backgrounds: flat 
spaces with zero-field or constant field. We use the zeta-function regularisation 
of determinants of second-order elliptic operators. The Dirac operator, If), is first 

order, but we shall define det Jp = J det Jp 1 . In order that Jp 1 may be an ellicp- 
tic operator, we shall work in Wick rotated Euclidean spacetime. In section 2 we 
consider (non-)Abelian gauge fields defined on multidimensional tori The 
non-trivial topology of these tori allows for the possibility of harmonic gauge po- 
tentials, that is, potentials which have vanishing field strength but which are not 
pure gauge. The existence of these harmonic gauge potentials is associated with 
the possibility of encountering non-trivial Wilson loops. On multidimensional tori, 
such gauge potentials are constant (up to a gauge transformation), thus allowing 
explicit construction of the eigenspectrum, zeta function, and effective action. The 
dependence of the effective action on these nontrivial Wilson loops may be viewed 
as a generalized Aharonov-Bohm effect. In sections 3 through 7, we consider gauge 
fields with constant field strength in arbitrarily many dimensions, working our way 
up from two dimensions. In all these cases we shall determine the eigenspectra of 
the gauged Laplacian and Dirac operator and shall calculate the associated zeta 
functions explicitly. Using special properties of zeta functions may often give the 
determinant and effective Lagrangian density in closed form. 

We shall also discuss the "folk theorem" C e g ~ B d / 2 \nB, and will point out 
a number of situations in which it is vio late d. The physical implications of the 
logarithmic term are discussed in reference H^l We shall show that this folk theorem 
is generally true in even numbers of dimensions, though there are exceptions, such 
as the scalar particle in two dimensions. In odd numbers of dimensions however, 
the logarithmic term is absent, and generically we obtain C e s ~ B d / 2 . 
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2. Field Free Gauge Potentials 



In this section we shall see that a generalized Aharonov-Bohm effect can influence 
the one-loop effective action of a system, even though the electric field strength is 
everywhere zero. Notice that an Abelian gauge potential is gauge-equivalent to the 
sum of its harmonic, monopole and coexact pieces 

A = H + A M + 5 X . (1) 

The number of independent harmonics equals the first Betti number and is thus 
intimately related to the topology of the spacetime manifold. We shall be interested 
in gauge potentials on d-dimensional tori, in which case there are d independent 
harmonics which may be chosen to be constant one-forms. In this section we shall 
further specialise to the zero-field case so that the monopole contribution and x 
are both zero. 

We begin with (constant) harmonic potentials defined on the d-dimcnsional 
torus R d /A, where the lattice A consists of points of the form Lj^^Ej, the 
n* being integers and the Ej being d linearly independent vectors in R d . It is, 
of course, possible to set L = 1 but we prefer to keep this prefactor in order to 
track the effect of scaling the lattice. Note that a constant field-free potential is 
not pure gauge since the "would be" gauge transformation is not single valued. 
These potentials are closed but not exact. Equivalently one sees that the Wilson 
loop W("fj) = exp(— i<f^ A ■ dx) = exp(— iA ■ Ej L) evaluated on the closed non 

contractible loop "fj from to LEj is gauge invariant, and is therefore an obstruction 
to gauging the potential away to zero. 

The eigenvalues of — D 2 and If) are in fact identical for zero field. On the 
multidimensional torus with periodic boundary conditions, ip(x) = ^(x + jl), ft £ A, 
the eigenvalues are easily computed to be 

A ™ VT/ ■ 9 lj (n t ~ ai)(n 3 - aj); n£Z d , (2) 



where the matrix (g 1 -*) is the inverse of gij — Ei ■ Ej] and di = ^-Ei-A. To compute 



L 

,J i ' -L^j , ctxiu Hi — 

the corresponding effective action one defines 

d 



5 cff = ilndet(- J D 2 ) = -i -(( 



(3) 



where (^(s) is the zeta function associated with — D 2 , 



C 



It should emphasised that the zeta function definition embodied in ^ implies 
both a regularization and a renormalization. When comparing the zeta function 
result with effective actions calculated using other renormalization prescriptions 
(e.g., Schwinger's) one should always bear in mind that different renormalization 
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prescriptions will yield effective actions that may differ by a finite renormalization. 
In particular, the quantity [i appearing in ((4|) is a renormalization scale, which has 
been introduced to keep the zeta function dimensionless, and the dependence of 
S c g on (i corresponds to a finite renormalization. The dependence on the renor- 
malization sc ale is logarithmic and proportional to the conformal anomaly (see 
e.g., reference ^^J) . In order that dimensionful quantities may be properly displayed 
throughout, we shall retain this normalization scale in all our calculations. 

In the present case the zeta function associated with the differential operator 
— D 2 is given in terms of a relatively well understood special function by 

£ (s)=2 (sr^ (s,5) - (s) 

The prefactor 2 arises from summing over both particle and antiparticle states, 
while the symbol (^ E (s,a) denotes the generalized Epstein zeta function, l -*-7 | 18 | 
defined by the sum 

Ci?( s >«) = ^I [g 13 (n% - ai)(rij - aj)] / (6) 
z d 

(the prime indicates that one should sum over non-zero eigenvalues only). 

The generalised Epstein zeta function is difficult to handle explicitly, at least in 
higher dimensions. Fortunately, we can apply the generalised Poisson resummation 
formula 

exp [— ng 10 ' (m — CLi){rij — aj)] = \j det gij ■ exp [— ngijm % vn? — 2irim l ai\ . 
z d z d 

(7) 

Taking a Mellin transform of the above gives | 1^ | 19 | 

C b m = ^r$r \/ det[ ^ ] E'Gfc" 1 ^')*"* e ~ 2 ™ mV (8) 

The zero-mode, = 0, is eliminated because the zeta function is defined by ana- 
lytic continuation in s. For large s the zero-mode makes no contribution. Equation 
ijHJ) is a generalisation of the well-known functional equation (reflection formula) 
for Riemann's (or rather Euler's) zeta function. ^Ql ft j s now relatively simple to 
show that 

( E (0,a) = (9) 
£ E (0,a) = T(d/2) Tr-^detfe] E'C^'^r 1 e^V 

The special case d — 2 has been extensively studied by Kronecker, ^] and for that 
reason we shall refer to „(0, a) as a generalised Kronecker sum. 
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Since E vanishes at the origin the prefactor (2tt/L) 2 
to the effective action and 

S cff = ±lndet(-.D 2 ) = -CsM- 



in (O does not contribute 



(10) 



Note in particular that in this case the effective action is independent of the renor- 
malization scale \x. This is actually a rather deep result, relate d t o the vanishing of 
the confbrmal anomaly for zero electromagnetic field strength. If we consider If) 2 
instead of —D 2 , the only change is to multiply the above with minus one- half the 
number of spinor components in d dimensions. This general result may be related 
to known results in one and two dimensions. 

In the simplest case, d= 1, the sum in (|9|) is a familiar trigonometric series. ^21 
We find the effective action 

LA" 



S eS = ilndet(-LH) = In 4 sin 



(11) 



Indeed, in one dimension one does not need the reflection formula. The exact zeta 
function is just a sum of Hurwitz zeta functions and the effective action can be 
computed directly. Using other techniques, we may show that the effect of a 
mass term is to change the determinant to 



Scff = ilndet(-L» 2 + m 2 ) 



ln<U 



LA 



sinh 



mL 



(12) 



2 J ' V 2 

Considering the next simplest case, we observe that in two dimensions any lattice 



may be re-scaled to satisfy 



1 Re(r) 
Re(r) rf 



(13) 



where r is a complex parameter known as the Teichmuller parameter. The Kro- 
necker sum becomes 

-2iri{ma\ -\-na-2 ) 



C(o> 



Im(r) 



E 



(14) 



7r ' — ' \m + rn\ 2 

This sum can be expressed in terms of Riemann theta-functions. I l^ l ^ l We find the 
effective action 



Cm 



= -2 log 



1 



t?(t) 



a i 

a-2 



(15) 



This two-dimensional 



Our theta-function conventions are those of Mumfo rd. 
effective action is in fact well known in string theory. SE1 

Though the discussion has, for clarity, been given in terms of an Abelian gauge 
potential, the extension to non-Abelian gauge potentials is simple. If the field 
strength is zero, we may use the non-Abelian version of Stoke's theorem to de- 
duce that the Wilson loops W(-f) — tr(Pexp § A ■ dx) form a representation of Hi, 

the first homology group. This implies that the A(x) = A yl a (:E) may be gauge fixed 
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to be mutually commuting constant matrices.^ The results of this section then con- 
tinue to hold provided one introduces an additional product over the gauge group, 

ndim(G) 
a=l 



3. Constant Field Strength: Two Dimensions 

Consider the case of a constant electromagnetic field in two dimensions. Let the 
area of spacetime be denoted by S. The field strength may be written as: 

'--(-bo)- <16 » 

As a potential we choose = (0, Bx) . 
Scalar Particles. 

With this choice the gauged Laplacian appropriate to scalar particles is 

- D 2 = -fl£- (d y -iBx) 2 . (17) 

By observing that — D 2 commutes with the momentum P y = —id y we see that on 
eigenstates of P v it reduces to 

-D 2 ^-d 2 x + B 2 ( X -^) 2 . (18) 

Since this is just (twice) the Hamiltonian of a harmonic oscillator it has eigenval- 
ues X Pt n = (2n + 1)|-B|. These eigenvalues are independent of p. Thus all levels 
are degenerate. Later we shall see that this degeneracy is just 2(\B\ • S)/2tv, the 
prefactor 2 again arising from the sum over both particle and antiparticle states. 
For notational simplicity take B to be positive. 

The zeta function of — D 2 + m is given in terms of a Hurwitz zeta function by: 



cm - ^ t r +i >f + "i ~* - ^ m - * + a- («> 



2vr \fi 2 J 



2ir 

n=0 

The one-loop effective Lagrangian density is £ c ff = <5'cff/S' = |lndet(— D 2 + 

m 2 )/S — —5^* {0)/S. Using properties of the Hurwitz zeta function discussed in 
the appendix, in particular its value and slope at the point s — 0, yields: 

£ cff = • m (^) - 1 - in (r ( i + ^)/V5F) . (20) 

It is easy to see that the m — > limit is well behaved — C e s — ► (_B/47r) • In 2 - 
while in the strong-field limit 

^ = -S-'°(7l) + s" l2 + 0<1) - (21) 

This simple example is already a counterexample to the folk theorem C c e ~ i? d / 2 In S. 
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The weak-field B — > limit may be taken by making use of the doubling formula 
for the Hurwitz zeta function: 



(22) 



which follows from the definition (|A.1[) of the appendix. Using this doubling formula 
we may write (|19[) as 



C(s) 



2BS 
2tt 

2BS 
2tt 



2B\ 
B 

7/2 



m 
' 2B 

771 



2B\ / 717r 

-T Ch [s; — 



2 

2B 



(23) 



With the help of equations (|A.9[) . (|A.11|) . and (|A.18|) , we may develop the weak- 
field expansion for the effective Lagrangian: 

[there is a typo, an extaraneous 2, in the second line of the published version] 



£ ff 



iC'(o) 


2 


S 


B 




2^ 





9 

m 
"•IT 



-In 

B 

"2n 

In 



(h[0 



2B 



-Ch 



2B 



to 

Ch(0;- 



Cff 0; 



2B 



B 



2k+2 



fc=l 



(2fc + 2)(2fc + l) 



to 
" ~B 

1 to 2 
2^ + 

2fc+l 



1 B 

12^2 



1 2B 
12^2 



7%) {2 



2fe+l 



-l} + 



B 



1 

2^ 



1 2 
— TO 

2 

n-1 



1 - In 



in 



B21C+2 



^ (2fc + 2)(2/c 



1 



" 12^2 

T, \ 2fe+2 



2fc+l 



- 1} +0 



B 



(24) 



Here the symbol B n denotes the nth Bernoulli number. Notice that there are no 
terms logarithmic in B in the weak-field expansion. However, as foreshadowed, 
the zeta-function renormalization has introduced nonstandard finite terms into the 
effective Lagrangian. An additional finite (often /i-dependent) renormalization is 
needed to remove these terms. Removing the constant term in the above equation 
corresponds to renormalizing the effective cosmological constant, while removing 
the term proportional to B 2 corresponds to a finite renormalization of electric 
charge. (We have chosen our notation in such a manner that electric charge does 
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not appear explicitly.) Having implemented these additional finite renormalizations 
we (finally) display the renormalized effective action as an asymptotic series starting 
with B 4 : 

(We have gone through this simple example in admittedly tedious detail because 
the same techniques will be used over and over in the following discussion.) 

In order to confirm the degeneracy factor we consider the heat kernel of — D 2 + 
m 2 . This heat kernel is found to be 

K(t) = tr' {exp([D 2 - m 2 ]i/^ 2 } 

OO 

= • e" 1 " 2 *^ 2 • cosech(Bi/M 2 )- (26) 
2tt 

Note that as t — > 0, K(t) — > 2(S'/i 2 /4nt), as it should according to the (known) short 
time behaviour of the heat kernel. This a posteriori proves that we have chosen the 
correct degeneracy factor 2(BS/2n). Equivalently, one may use the fact that as 
B — > the heat kernel must approach that of the free two-dimensional diffusion 
operator, to obtain an alternative verification of the degeneracy factor. It should 
be emphasized that we are considering the heat kernel appropriate to a complex 
scalar field. That is to say, the sum in equation ([26]) includes both particle and anti 
particle sectors. 
Dirac Particles. 

The analysis for Dirac spinors closely parallels that of the charged scalar field. The 
square of the Dirac operator is 

p 2 = -D 2 + Z^F^ = -d 2 - (d y - iBxf + l5 B. (27) 

The eigenvalues are A n = (2n + 1)\B\ ± B. We now re-index the states so that the 
eigenvalues become A„ = 2n|£?|, where the effective degeneracy is (\B\S/2tt) for 
n = and 2(\B\S/2n) for n > 0. 

The zeta function of Ip 2 is given in terms of Hurwitz zeta functions (henceforth 
we take B > 0), 

c w ^{,g(^p(-)-'} 
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The one-loop effective Lagrangian density for spinors is £ c g = +5^* (0)/S, so that: 

B + rn 2 [2B\ B , , 2, B , /27rm 2 \ , . 

l-(^j+^^(l + ^)--ln(-^j. (29) 

The weak-field limit B — > for spinors is obtained just as it was for scalars. As 
in the scalar case we must implement an additional finite renormalization of the 
cosmological constant and the electric charge. Then we obtain 



The limit m — > of equation (T2"9")l is ill behaved (£ e ff — > 00). This is an infrared 
singularity associated with the fact that the Dirac operator develops a zero-mode 
as m — > 0. One may deal with the zero-mode by simply dropping the ground state 
(n = 0) from the sum prior to taking the m — > limit. Equivalently, for masslcss 
Dirac spinors one must exclude the zero mode "by hand" . Recall that the definition 
of the zeta function is in terms of ^ not J^. The difference has up to now been 
irrelevant. Taking careful note of this difference leads to: 

f~ 1 \ 2BS^(2nB\ s 2BS f2B\~ s /■ 
The one-loop effective Lagrangian density is simplified to 

2ir V 71 "/ 1 



' — Bok+2 ( 2B 



C cS = +— ■ In ( — ) . (32) 



The heat kernel of Ip 2 + m 2 is readily evaluated 

, BS , ( Bt\ ( mH\ 
i^) = -coth(^Jexp(-^J. (33) 

Note that as t -> that if(t) -> 2 • (SfP/iirt), as indeed it should. (The factor of 2 
reflects the existence of two spinor components in two dimensions) . In the massless 
case the "zero-mode-suppressed" heat kernel is 



colli 1^)1 



(34) 



Note that K '(t) — > as t — ► 00, thanks to the explicit exclusion of the zero mode. 

Finally, note that calculation of the heat kernels allows one to deduce all the 
Seeley-deWitt coefficients for the case of constant field. Use of the Taylor series for 
cosech(Bt) and coth(Bt) yields: 

a 2w (-D 2 ) = -2 {Z l) B 2n -£ 2 ", (35) 

a2n(p 2 ) = [ -j^-B 2n -B 2n , (36) 
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where B n denotes the nth Bernoulli number. The 'odd' coefficients a2 n +i all vanish. 
It is very unusual to know all the a n . For arbitrary fields these coefficients are leading 
order terms in a low-momentum approximation. Thus they remain interesting for 
general fields. 



4. Constant Field Strength: Three Dimensions 

In three dimensions the field strength may be written as: 




F,j,v = I B (37) 
\0 —B / 

Let the volume of spacetime be denoted by V = S ■ L. In the gauge = (0, 0, By), 
we may immediately write down the three-dimensional heat kernel Kd=z{t) 
Kd=2(t). Three-dimensional zeta functions quickly follow: 

r(s) J 



/47Tt 



r(«) 



■C= 2 (*-|). (38) 

All factors are analytic at s — 0. In particular, since j^jy ~ s + 0(s 2 ), the derivative 
at s = is 

Ca(°) - ^ ' r (-i) ' C= 2 H) = ■ C= 2 H)- (39) 

This means that the one-loop effective Lagrangian densities are given in terms of 
^ „(- |,x). Using results of the two-dimensional discussion we see: 



Scalar Particles: £ c ff = • (2B) 3/2 . £ H (-*;*+ 

Dirac Particles: £ cff = -— • 2B ■ ^zV2B C H (~b 1 + Sf) + m } ■ ( 40 ) 



8tt 

Note that these effective actions are independent of the renormalization scale /i. This 
(nonobvious) result is a consequence of the vanishing of the conformal anomaly in 
odd-dimensional space-times. For strong fields, the Hurwitz zeta function can be 
computed by convergent series. [See equation (|A.6|) of the appendix.] For the cases 
of interest: 

C (-!• i + ^) = C (-i^)-T(-) l (2l ~ m -(—\ -C r-i + z-i) 

SffV 2>2^2B> SifV 2^2) Z^ V ' O l V \ 2B ^ H 2 Tt '2<'' 

1=1 • V / 

f=l ■ V / 

The coefficients in these expansions can be obtained (with severely limited accuracy) 
from the tables of Jahnke and Emde. ^1 Alternatively, the proliferation of personal 
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computers allows these coefficients to easily be computed to any desired accuracy. 
Note the presence of fractional powers of the field strength and the absence of 
logarithmic terms when equation (|41|) is substituted into (|40j) in order to obtain 
an expansion for the effective Lagrangian. We shall see that these general features 
persist in any odd number of dimensions, at least for constant field strength. The 
leading strong-field limit in three dimensions is 

Scalar particles: C cff = — (2B) 3/2 ( H 



Spinor particles: 



£cff 



47T 



I. P 

2' 2) 



(42) 



The leading term in the strong-field expansion for the spinor Lagrangian has pre- 
viously been obtained by Redlich. 1^3 

The leading coefficients appearing in the strong-field limit have been evaluated 
numerically. This was done by using the reflection formula for the Riemann zeta 
function to write C, H {— \\ 1) = C r(~ I) = — 5tCr(§) - Note that R (z>) is given 
by a nicely convergent series suitable for computer evaluation. We find R {\ ) « 
2.612375, Cfl(~5) ~ -0-207886. In a similar vein, we use the "doubling formula", 
( H (s: i) = (2 s - l)( R (s) to deduce 



l. i* 

" 2 ' 2 ' 



1 

47T 



' -£}C„CI> 



0.060888. 



(43) 



In order to obtain the weak-field limit for the effective Lagrangian, we return to 
equation (f40|) . We shall need equation (1A.7|) of the appendix, evaluated at s = |. 



'rV r(-i) 



x 3 / 2 L 

r(2fc 



2 V 2 



2/* 3/2-2fc 



fc=l 



(2fc)! 



■O(o: 



l/2-2n 



(44) 



After an additional finite renormalization of the cosmological constant and the 
electric charge, the weak-field limit for the scalar effective Lagrangian is 

2k+2 



c 



efi 



m 

7T 



r(2fc + i) 



2k+2" 



'(2fc + 2)!P(-± 
while the spinor Lagrangian is 

7 Hoi- i o 

4vr ^ 

fe=l 



>-''■• 



2k+2 



(2fe + 2)!r(-|) 



25 

2 



2fc+2 



2n+4^ 



0(B 



2ri+4\ 



(45) 



(46) 



Only even, positive-integer powers of the field enter the weak-field expansion. 

Notice the absence form our discussion of the Chern-Simons secondary char- 
acteristic class (the topological mass term for the photon). This term might be 
expected to appear on rather general grounds, following arguments of Niemi and 
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Semenoff ESI anc j Redlich^ll. However, a simple analysis is s ufficie nt to show that 
the Chern-Simons term vanishes for constant field strength. 1 27 1 29 1 w e should also 
note the appearance of fractional powers of the field strength has also been noted 
in Redlich's calculation for a massive spinor. L ^The present calculation generalizes 
this result to the massive case and also to scalar particles. 



5. Constant Field Strength: Four Dimensions 

In four dimensions the field strength may be written as: 



F„ 



/ E 0\ 
-E 
B 
V -B J 



(47) 



so that the four-dimensional problem breaks up into two two-dimensional ones. 

This result follows from the block-diagonalizability of antisymmetic matrices 
through orthogonal transformations. It is important to know that E and B are 
invariant scalars that characterise the electromagnetic field. Note that E 2 + B 2 = 
\F^ V F^ V = \F 2 , and 2EB = e w F"" 'F ap = FF. While it is more common 
to characterise the field in terms of F 2 and FF, it is more useful for us to use 
the invariants E and B. Though we shall be working in Euclidean space, we note 
that the result (|4"T|) also obtains in Minkowski space. This follows from the fact 
that it is always possible to make a Lorentz boost such the electric and magnetic 
fields become parallel and then (by a rotation) to make both fields point in the x 
direction. 



5.1. E — 0, B ^ 

Our analysis in this case parallels that of three dimensions. Writing the volume of 
spacetime as ft = ft± ■ S, and choosing the gauge as = (0, 0, 0, By), we factorize 
the heat kernel 

K d ^{t) = D^^- ■ K d=2 {t). (48) 

47TE 

Here D = 1 for scalars and D = 2 for spinors, reflecting the fact that four- 
dimensional spinors possess twice as many degrees of freedom as two-dimensional 
spinors. The four-dimensional zeta function is given by 

Inserting the results of the two-dimensional case yields 

B 2 n (2B\~ S (, H {s-l 



Scalar Particles: £(s 



' 2 ~ IB > 



2tt 2 \fi 2 J (s - 1) 



June 16, 2009 6:34 WSPC/Guidelines-IJMPA analytic-4-arXiv 



Analytic results for the effective action 13 



Dirac Particles: C (s) = —s- • — • ' 2B 

3 TT J \ fl d ) \S — 1) 



m 2 BQ ( (m/fi) 



4ir 2 1 s - 1 
Taking derivatives at s = is straightforward. 



(50) 



5.1.1. Scalar Particles: 



For the effective Lagrangian density one obtains 
B 2 



4ir 2 



, 2B 

1 - hl 1 7? 



O-ijs + SJ + ChHjs + w)} 



(51) 



Using results from the appendix one finds that 



B 2 

24 



B A 



■IB 



\n{T{\+y)/V^)dy. 



(52) 



This form makes it easy to extract the large field limit B ^> m 2 , 



47T 2 



B 2 
24 



m 
~8 



Bm 2 



— — In 2 



(53) 



The coefficient H (— 1; ^) has been evaluated numerically. This was done by noting 
that H (s; i) = (2 s — 1) • ^ R (s); and by using the reflection formula (see appendix) 

to relate £ H (-l; 5) to C^(+ 2 )- We obtain Ch^ -1 ' i) ~ °- 053829 - 

In order to obtain the weak-field limit it is easiest to return to equation (|51ll 
and to use the doubling formula (]22[) along with equations (|A.10j) and (|A.16|1 of the 
appendix. This procedure yields an asymptotic expansion in terms of powers of B 2 ; 
in particular there are no terms logarithmic in field strength. After performing the 
usual additional finite renormalization of cosmological constant and electric charge, 
we obtain the extended Euler-Heisenberg effective Lagrangian 



B 2 
16tt 2 



-7B 2 
360^4 



E 



B 



2k+2 



< (2fc + 2)(2fc + l)(2fc) 



[2 2k + 2 - 2] 



B 



2k 



0(B 



2n+2\ 



(54) 



This result should be compared with the Minkowski-space Euler-Heisenberg La- 
grangian obtained by Schwinger.^In order to make the comparison, we must take 
Schwinger's formula and after setting E = multiply by — 1. In the general case, for 
which E does not vanish, we must replace the £J's in Schwinger's formula with iE 
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before multiplying by —1 to obtain the Euclidean result. Finally we have implicitly 
chosen units for the renormalized charge so that Schwinger's a 2 corresponds to our 
1/167T 2 . After making these adjustments for notation and metric signature, we may 
confirm that our expansion agrees with Schwinger to the leading order (which is 
the only order explicitly displayed in reference Note though, that in this zeta 
function formalism, it is easy to display all orders of the asymptotic expansion for 
the Euler-Heisenberg Lagrangian in the weak-field limit. 



5.1.2. Dirac Particles: 



For Dirac spinors one obtains 



C 



off 



B 2 
"2tP 



1 - * (f )] i + Si) + C(-i; i + £)} 



+■ 



n 2 B r 
8tt 2 [ 



hi 



(55) 



A little work using results given in the appendix yields 



2tt 2 



12 



m 2 B 



m B 



1 - In 



2B 

,,2 



5 2 C(-i) 



m 2 B 





- 1 


-B 2 [ 






Jo 



2B 



ln(r(l + y)/V2^)dy 



(56) 



For strong fields one uses lnT(l + e) = —76 + 0(e 2 ), to establish 



Ce 



1 

2^2 



B^_ 
12 

m 2 B 



m 2 B 



. m , 



1 - In 
m 2 B 



2B 

,,2 



ni 
~~8 



m 2 B 



| ln(27r) + |m 4 \ + m A 0(m 2 /B). (57) 



The coefficient ^ R (— 1) ~ —0.165421 has been calculated numerically. 

In order to obtain the weak-field limit we return to equation (|55[) and employ 
the identity Ch(s', 1 + x) = Of ( s ! — ^~ s - Just as in the scalar case we obtain an 
(extended) Euler-Heisenberg Lagrangian. 

2k 



■Ceff = 



B 2 
16^2 



45 m 4+ ^ 



4 B 



2k+2 



^ (2fc + 2)(2fe + l)(2fe) \m 



(58) 



which specializes to the Euclidean version of Schwinger's result. ^ 

The results of this section may easily be extended to discuss the case B = 0, 
E 7^ 0. One needs merely replace J3 by E in all formulae to obtain valid Euclidean 
space results. Recall, though, that in continuing to Minkowski space one effects 
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the transformation E — > iE. In the next section we shall begin to address the 
complexities encountered when E and B are both nonzero. 



5.2. E y£ 0,B ^ 

The scalar spectrum is given by A = (2n + 1)\E\ + (2n' + 1)\B\ + to 2 , where all the 
states have the degeneracy 2{\E B\£l) / At: 2 . The zeta function is 



\EB\n 

2tt 2 



E 

n,n/=0 



(2n+l)|£| + (2n' + l)|B|+TO 2 



(59) 



For spinors the analogous zeta function is a trifle more complicated [the published 
version of this equation has a typo]: 



\EB\n 
~^ 2 ~ 



4- E 

n,n' = 1 
oo 

+2-E 

n=l 

oc 

+2-E 



n'=l 



2n|£| +2n'|S| 



2n\E\+m 2 

2n'\B\ + to 2 
M 2 



+1- 



r 2 - 




m 









(60) 



These zeta functions are in general so unwieldy as to be unmanageable. To proceed 
we restrict ourselves to the case \E\ = \B\ = F, corresponding to self-dual and anti- 
self dual constant fields. In this special situation we can recast these zeta functions 
in terms of the elementary Hurwitz zeta functions. 
For the case of a scalar particle 



Cm 



F 2 n 

2tt 2 



(Tf - E {n + n' + l + ^y S 

V ^ ' n,n'=0 



(61) 



We define t = n + n', and observe that there are precisely t + 1 ways in which n 
and n' can be arranged to sum to t, consequently 



C(«) 



F 2 n 

~2^ 



F 2 n (2F 

~2^ ' \^ 

F 2 n 



2tt 2 



(9 TP\ ~ s 00 — 

=JJ •E^ + 1 )-(*+ 1 + ^) _ " (62) 

{Q S -l;l + ^ ) -^.Q S ;l + ^)} (63) 

{C ff ( s -i;fF)-iF-C ff ( s ;fF)}- (64) 
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Here we have used (h(s; l + x) = (h(s; x) — x s . For the Dirac spinor, very similar 
results obtain, though one must be careful with degeneracy factors: 



Coo 



F 2 n 

4tt 2 ' 



2F 
2F 



4 £ (n + n'+^Y 



Re-indexing the sum, using t = n + n' — 1, yields 

F 2 il /to 2 



Cs P inor( S ) ~ 2 CscalarM 



47T 2 \ fl 2 



(65) 
(66) 

(67) 



Apart from the ground-state contribution (n = n' = 0; A = to 2 ) the scalar and 
spinor zeta functions are proportional to one another. This desirable property does 
not, unfortunately, continue to hold if \E\ ^ \B\. To discuss the effective action, the 
strong-field limit, and the weak-field limit, it suffices to discuss the scalars. 
The effective Lagrangian for scalars is 



2F ^5 H 



Using results from the appendix yields 







'2F\ 


—F 2 - -m 4 


In 




12 8 







4tt 2 

i^ 2 
1 

i^ 2 ^ 

^■JJ* ln(T(l+y)/V^)dy. 



r(i + w) 



2^ 



(68) 
(69) 

(70) 
(71) 
(72) 

(73) 



The strong-field and weak-field limits may be written down in the standard manner, 
we omit details and present the results. For strong fields 



i 



+ \m z F + -(1 + 7 )to 4 + O ( — 



(74) 
(75) 
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For weak fields, the Euclidean-space (scalar) Euler-Heisenberg Lagrangian is 



16tt 2 



B 



2k+2 



-IF 2 " 

15 m4 + hi ( 2fc + 2 )(2/c + l) 



4 + ?U?? 



2/,- 



OF 



2rc+2\ 



(76) 



When continued to Minkowski space, the lowest term in this expansion reproduces 
Schwinger's result. ^ To relate this scalar result to the case of a massive Dirac 
spinor, one notes that 

^•ln(m 2 /^ 2 ). (77) 



C 



Spinor 



= -2-L 



Scalar 



The massless case affords considerable simplifications; for the scalar one obtains 

F 2 fl (2F S 



C(«) 



2tt 2 



This leads to an effective action 

"47T 2 ' 



■Ccff 



15 ln 



1) 



2F 
1? 



(78) 



(79) 



The term proportional to F 2 is an artifact of the zeta-function method, and may, 
in the usual fashion, be removed through a finite renormalization of the electric 
charge. For massless Dirac spinors, explicit exclusion of the zero-mode leads to 
^Spinor = — 2£scaiar- These massless (anti-)self-dual effective actions have been 
previously discussed in the literature, see for example | 15 | 30 | 31 ] 



6. Even Dimensionality 

In d = 2N dimensions a constant field may be brought into the block diagonal form 

( B x \ 
-Bx 



F, 



V 



(80) 



B N 
-B N / 

where the ±Bj are the zeros of the characteristic polynomial det(B ■ I + iF). The 
problem thus decomposes into N two-dimensional problems, allowing the associated 
zeta functions to be written down by inspection. 
Scalar Particles 

m 2 with degeneracy d n = 



{Ef=i(2^ 

Ilili 2(|^| • Si/2n). The zeta function is 



The eigenspectrum is 



1)1-8,1} 



Co 




m 2 



m 



(81) 
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-^-r-{n(^)«-(^))„ P( -f,. m 

Dirac Particles 

The eigenspectrum is = {J2iLi ^ n i\Bi\} + m 2 with degeneracy d n — YliLi(\^i\ ' 
Si/2n) ■ 2 a '(™*),where w(0) = 0, and ui{n > 0) = 1. The zcta function is 

Cw-n(^)-E^'( {£L2 ";'f" + "'T ' <»> 





«»M^I >^{-^p) < 84 > 

These zeta functions are tedious to calculate with in general, to proceed we make 
the radically simplifying assumptions that \B\\ = \Bz\ = . . . = B, and that m = 0. 
Of course, in two dimensions the condition on the field is always met since there is 
only one B. In four dimensions it means that we confine ourselves to (anti-)self- 
dual backgrounds. With these assumptions the scalar zeta function is simplified 
to 

^■©"■T'f^- hf) " (85) 

while for the massless spinor, exclusion of the zero mode leads to 

cw-(^)"-^-r^««* t »"'- 1 )- <«> 

Note that for N = 2, (i.e., four dimensions), the spinorial zeta function is four times 
the scalar zeta function, in agreement with our earlier result. 

For these backgrounds we can relate the zeta function in d-dimensions to the 
ones in (d — 2) and (d — 4) dimensions by recursion relations. For the scalar case, 
a repeated integration by parts establishes 

so that for the scalar there are two disconnected series starting with d = 2, d = 4, 
respectively. For the spinorial case a single integration by parts suffices to obtain 



WSj ~ 2tt ' N-l \27r ' Sw - 2 



(a). (88) 



In this case there is a single recursive series. We define N=0 (s) = 0, while our 
previous calculations have shown N=1 (s) = ^p- • (2B/fi 2 )~ s ■ ^ R (s)- So we de- 
duce C, N=2 ( s ) = ("ir) 2 ' • (^ R (s — 1), which verifies our previous four 



' N=2 

dimensional calculation. 
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Without resorting to numerical computation, we can make some general ob- 
servations regarding the d-dimensional case. Note that all the (i-dimensional zeta 
functions are of the form 

CwatW = Bd/2 ■ (^) (89) 

in which case 

C cS ex B d / 2 {ln(B/(x 2 ) ■ /(0) - /'(0)} (90) 

so that the "folk theorem" generically holds for this class of fields. The "folk the- 
orem" can fail if /(0) happens to be zero; this is in fact exactly what happens for 
massless scalar particles in two dimensions. We remind the reader that the integrals 
flSSJ) and (j5B"|) defining f(s) make sense only for Re(s) > d/2. The value of /(0) is 
defined by analytic continuation in s. 

We also wish to point out that our recursive formulas (f87|) and (|88|1 are not the 
only way of proceeding. The spectra are in this case sufficiently simple that explicit 
calculation in terms of Hurwitz zeta functions is possible. 



7. Odd Dimensionality 

In d = 2N + 1 dimensions, we write the volume of spacetime as fid — L ■ Q,2N- The 
heat kernel is related to that in 2N dimensions by Kd{t) — -7= ■ K2N(t), so that 
(paralleling the discussion of the three-dimensional case) one deduces 

0«) = ^-^'C(»-|). CD 

Consequently £ d (0) = —L ■ £ 2 jv( — f)' an< ^ * n ^ e nota ti° n of equations ((89| and 
(|90| one has 

C cS o,B d ' 2 ■ f{-\). (92) 

Though this calculation is carried out for the special case \B\\ = ■ ■ ■ = \Bi\ = ■ ■ ■ = 
I -Bw I = B, we expect that for general fields the result C c s ~ B d l 2 will remain true 
in any odd number of dimensions. 



8. Conclusions 

We have studied the effect of constant field strengths on the determinants of the 
Dirac operator and gauged Laplacian. We began with the topologically interest- 
ing cases of harmonic gauge potentials, before proceeding to monopole potentials 
with constant field strengths. In both these cases we obtained the spectrum of the 
squared Dirac operator and gauged Laplacian explicitly and calculated the corre- 
sponding zeta function directly. Computation of the effective action then reduces 
to a relatively simple application of special functions. 
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Appendix A. Hurwitz zeta Functions. 

To make this paper as self contained as possible, we shall present here a number of 
definitions and useful results related to zeta-function theory. Most of these results 

IOQI 

may be obtained by consulting Gradshteyn and Ryzhik. ° Other useful references 
include Weil^^, Abramowitz and StegunJ^The Encyclopedia of Mathematics, ^31 
and the Bateman Manuscript Project. <^ 
The Hurwitz zeta function is defined by: 

oo 
n=0 

This series converges absolutely for Re(s) > 1, and the function so defined may be 
analytically continued to the entire complex plane. There is a single simple pole at 
s = 1 and in that neighborhood H (s, x) = j^j — ip(x) + o(s — 1). Setting x = 1 
reproduces the ordinary Riemann zeta function, 

oo 

C*00 O-i: X> ( A - 2 ) 

71=1 

The Riemann eta function is closely related to the Riemann zeta function. It is 
defined by 

oo 

^) = E(~ 1 )" +1 "^- (A-3) 

71=1 

The Riemann eta function, being defined by an alternating series, is numerically 
much better behaved than the Riemann zeta function, and furthermore converges 
over a larger region of the complex s-plane [Re(s) > 0]. In terms of the eta function 

C^HttV^ (s) - (a - 4) 

The Riemann zeta function satisfies the classical reflection formula 

O*) = ^ • My ) ' r(l - S) ■ ( R (l S). (A.5) 

By combining the reflection formula with an improved series representation for 
the zeta function given in terms of the eta function one may numerically compute 
Riemann's zeta function over the entire complex plane. 

The Hurwitz zeta function may be evaluated in terms of Riemann's zeta function 
by using the binomial series to write 



-s-i i 
n ■ x 



( H ( S ;i + x) = J2J2 

71=1 1=0 

= ( R (s) + t(-) 1 • * > + 1 >--;/' + Z - 1 > • ( R (s + l)- X K (A.6) 
1=1 



June 16, 2009 6:34 WSPC/Guidelines-IJMPA analytic-4-arXiv 



Analytic results for the effective action 21 



This series is con vergent for |x| < 1. An asymptotic expansion for large x is provided 
by reference 1251 [equation 1.18 (9)]: 



1 



r( s ) 

n 

fc=i 



^- s r(s - 1) + 5x" s r(s) 

T(s + 2k + l) ,_ 2 _ 2fc _ X _ S _2„ 



2k" 



(2fc)! 



0(af 



(A.7) 



At special values of s more information is available. The value of the Hurwitz 
zeta function at non-positive integers is known in terms of the Bernoulli polynomials 

B n+1 (» B^ +2 (a;) 



n ' ,X) n + I [n + [)(,, -!-2)' 

in particular, 

H (— 1; x) = — + 2 X ~ 12 = — 2 C 35 — 2) ^~ 24' 

The derivative of the Hurwitz zeta function at s = is known to be 
£ H (0;x) = \n(T(x)/V2^). 



(A.8) 

(A.9) 
(A.10) 

(A.ll) 



Derivatives at other values of s are not given in the standard tables. One may make 
some progress by noting that 



Ox 



■( H (s + l;l + x) 



(A.12) 



so that 



(A.13) 



d( H (s;l + x) = _Q g +1 . 1 + x) _ s .£ h{s + 1]1 + x) 
This recursion relation, when applied to the case s = — 1 yields 

C(-!; 1 + x) = Cui' 1 ) + H* 2 + x) + [ ln(r(l + y)/V^)dy 

Jo 

= Cxi- 1 ) + W 2 + x ) - 3-^ln(27r) + xlnr(l + 

yip(l + y)dy, (A.14) 



This is the best analytic result that we have been able to obtain. For |x| < 1 we 
may make the convergent expansion 



C' ff (-l;l + .x) = O- 1 ) + W 2 + x ) - \x\xi{2tt) + xhxT{l + x) 



k=2 



(fc + 1) 



(A.15) 



June 16, 2009 6:34 WSPC/Guidelines-IJMPA analytic-4-arXiv 



22 S. K. Blau, M. Visser & A. Wipf 



For |x| large, a computationally usef ul asymptotic expansion may be found in ex- 
positions by Elizalde and Soto^SEIl 



*y £[( 1'^) (2*^ 2 ^ 12) ' ^ n ^ 4*^ 12 



ri-1 td 

(But note theat the final sign is erroneously displayed in reference ^l) We have 
resorted to numerical methods to evaluate 

f a (§) pb 2.612375; Cb( - s) ~ -0-207886; 
(^(2) w -0.937548; Cj( _1 ) ~ -0.165421. (A.17) 
Finally, we note that r(^) = y/ir, and that Stirling's approximation is 

m(r(,)/^) = ^)- +l2 -+g (2fc+2 2 ; ( ;: +1) ^ i - 2fc + o(,- ) . (A.18) 

It is often sufficient to use the simpler form 

ln(r(a;)/v^7r) ={x-\) ln(x - 1) - (a; - 1) + 0(l/x). (A.19) 
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